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Individual metallic single-wall carbon nanotubes show unsual non-Ohmic transport behaviors at 
high bias fields. For low resistance contact samples, the differential conductance dl/dV increases 
with increasing bias, reaching a maximum at ~ 100 mV. As the bias increases further, all /dV drops 
dramatically [Yao et at, Phys. Rev. Lett. 84, 2941, 2000]. The higher the bias, the system behaves 
in a more normal (Ohmic) manner. This so-called zero-bias anomaly is temperature-dependent (50- 
150 K). We propose a new interpretation. Supercurrent runs in the graphene wall below ~ 150 K. 
The normal conduction-electron currents run outside the wall, which are subject to the scattering 
by phonons and impurities. The currents along the tube induce circulating magnetic fields and 
eventually destroy the supercurrent in the wall at high enough bias, and restore the Ohmic behavior. 
If the prevalent ballistic electron model is adopted, then the scattering effects cannot be discussed. 

PACS numbers: 73.23.-b, 73.63.-b, 72.80.Vp, 74.78.Na 



I. INTRODUCTION 

In 2000, Yao, Kane and Dekker [l[ reported the 
high-field transports in metallic Single- Wall carbon Nan- 
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FIG. 1: Typical current I and differential conductance dl/dV 
vs voltage V obtained using (a) low-resistance contacts (LRC) 
and (b) high-resistance contacts (HRC). The inset to (b) plots 
dl/dV vs V on a double-log scale for the HRC sample. After 
Yao et al. [l|. 



oTubes (SWNT). In Fig. Q] we reproduced their I-V 
curves, after Ref. Fig. 1. At low fields (voltage 
~ 30 mV) , the currents show temperature-dependent dips 
near the origin, exhibiting non-Ohmic behaviors while at 
the high fields (~ 5 V) the resistance R versus the bias 
voltage V shows a relation: 



R = R Q + V/I (high VO, 



(1) 
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where Rq and Iq are constants. The original authors dis- 
cussed the low-field behavior in terms of one-dimensional 
(ID) Luttinger liquid (LL) model. Many experiments 
however indicate that the electrical transports in SWNT 
have a two-dimensional (2D) character. Q In fact, the 
conductivity in individual nanotubes depends on the cir- 
cumference and the pitch characterizing a space-curve 
(2D). Hence the nanotube physics requires a 2D theory. 
In the present work, we present a unified microscopic the- 
ory of both low- and high-field conductivities. Carbon 
nanotubes are discovered by Iijima. [3j The important 
questions are how the electrons or other charged particles 
traverse the nanotubes and whether these particles are 
scattered by impurities and phonons or not. To answer 
these questions, we need the electron energy band struc- 
tures. Wigncr and Seitz (WS) Q developed the WS cell 
model to study the ground state of a metal. Starting with 
a given lattice, they obtain a Brillouin zone in the fc-space 
and construct a Fermi surface. This method has been 
successful for cubic crystals including the face-centered 
cubic (fee), the body-centered cubic (bec), diamond and 
zincblende lattices. If we apply the WS cell model to 
graphene, we then obtain a gapless semiconductor, which 
is not experimentally observed. [2] We will overcome this 
difficulty by taking a different route in Sec. |TTJ 

SWNTs can be produced by rolling graphene sheets 
into circular cylinders of about one nanometer (nm) in 
diameter and microns (fJxn) in length. [5| The electri- 
cal conduction in SWNTs depends on the circumference 
and pitch, and can be classified into two groups: either 



2 



0.5 
T= 0.0 



-0.5 




-4 >2 c a A 

Biag voltage (mVJ 



FIG. 2: Current-voltage curves of metallic SWNT at gate 
voltages of 88.2 mV (trace A), 104.1 mV (trace B) and 
120.0 mV (trace C). After Tans et al. in Ref. [g, Fig. 2. 

semiconducting or metallic. [2] In our previous work Q, 
we have shown that this division in two groups arises as 
follows. A SWNT is likely to have an integral number 
of carbon hexagons around the circumference. If each 
pitch contains an integral number of hexagons, then the 
system is periodic along the tube axis, and "holes" (not 
"electrons") can move along the tube. Such a system is 
semiconducting and its electrical conductivity increases 
with the temperature, and is characterized by an activa- 
tion energy e 3 . Q The energy e 3 has a distribution since 
both pitch and circumference have distributions. The 
pitch angle is not controlled in the fabrication processes. 
There are far more numerous cases where the pitch con- 
tains an irrational number of hexagons. In these cases, 
the system shows a metallic behavior experimentally ob- 
served. [8fl 

We primarily deal with the metallic SWNTs in the 
present work. Before dealing with high-field transports, 
we briefly discuss the low-field transports. Tans et al. 
[i| measured the electrical currents in metallic SWNTs 
under bias and gate voltages. Their data from Ref. |9|, 
Fig. 2, are reproduced in Fig. [21 The currents versus the 
bias voltage are plotted in Fig.[2]at three gate voltages: A 
(88.2 mV), B (104.1 mV), C (120.0 mV). Significant fea- 
tures are: 

(i) A non-Ohmic behavior is observed for all, that is, 
the currents are not proportional to the bias volt- 
age except for high bias. The gate voltage charges 
the tube. The Coulomb (charging) energy of the 
system having charge Q is represented by 

Seoul - Q 2 /2C, (2) 
where C is the total capacitance of the tube. 

(ii) The current near the origin appears to be constant 
for different gate voltages V^ a te, (A)-(C). This fea- 
ture was confirmed by later experiments, [8lll0j The 



current does not change for small varying gate volt- 
age in a metallic SWNT (while the current (mag- 
nitude) decreases in a semiconducting SWNT). 

(iii) The current at gate voltage V ga .te = 88.2 mV (A) 
reverts to the normal resistive behavior after pass- 
ing the critical bias voltages on both (positive and 
negative) sides. Similar behaviors are observed for 
V gate = 104.1 mV (B) and V gatc = 120.0 mV (C). 

(iv) The flat current is destroyed for higher bias voltages 
(magnitudes). The critical bias voltage becomes 
smaller for higher gate voltages. 

(v) There is a restricted Vg a t c -range (view window) in 
which the horizontal stretch can be observed. 

Tan et al. [9| interpreted the flat currents near Vbias = 
in Fig. [2] in terms of a ballistic electron model [2j. 

We propose a differing interpretation of the data in 
Fig. [5] based on the Cooper pair [ll| (pairon) carrier 
model. Pairons move as bosons, and hence they are pro- 
duced with no activation energy factor. All features (i)- 
(v) can be explained simply with the assumption that 
the nanotube wall is in the superconducting state as ex- 
plained below. 

The supercurrent runs without obeying Ohm's law. 
This explains the feature (i). The supercurrents can 
run with no resistance due to the phonon and impu- 
rity scattering and with no bias voltage. Bachtold et al. 
[12j observed by scanned probe microscopy that the cur- 
rents run with no voltage change along the tube in 
metallic SWNTs. The system is in a superconducting 
ground state, whose energy —E g is negative relative to 
the ground-state energy of the Fermi liquid (electron) 
state. If the total energy E of the system is less than the 
condensation energy E g : 

E = K + E Coul + E^<E g , (3) 

where K is the kinetic energy of the conduction electrons 
and the pairons, and 

E^ = QV hias (4) 

is the Coulomb field energy, then the system is stable. 
Experiments in Fig. [2] were done at 5 mK. Hence, we may 
drop the kinetic energy K hereafter. The superconduct- 
ing state is maintained and the currents run unchanged 
if the bias voltage Vbias is not too large so that the in- 
equality ([3]) holds. This explains the horizontal stretch 
feature (ii). 

If the bias voltage is high enough so that the inequal- 
ity symbol in Eq. ([3]) is reversed, then normal currents 
revert and exhibit the Ohmic behavior, which explains 
the feature (iii). 

The feature (iv) can be explained as follows. For higher 
V^ a tc there are more amounts of charge, and hence the 
charges Qa, Qb, Qc for the three cases (A, B, C) satisfy 
the inequalities: 

Qa < Qb < Qc- (5) 



The horizontal stretches are longer for smaller bias volt- 
ages. At the end of the stretch (Vbi as ,max) the system 
energy equals the condensation energy E g . Hence, we 
obtain from Eq. ([3]) after dropping the kinetic energy K 

^,max — Q^biaSjinax = Q^max 

= E g - E Coul =E g - Q 2 /2C. (6) 
Using Eq. (JSJ, we then obtain 

Knax,A > Vmax.B > V maXt C, (7) 

which explains the feature (iv). 

The horizontal stretch becomes shorter as the gate 
voltage legate is raised; it vanishes when V ga ± c is a little 
over 120.0 mV. The limit is given by 

i^.iuax = Eg — -Ecoui = Eg — Q 2 /2C = 0. (8) 

If the charging energy -Ecoul exceeds the condensation 
energy E g , then there are no more supercurrents, which 
explains the feature (v). Clearly the important physical 
property in our pairon model is the condensation energy 

Eg- 

In the currently prevailing theory [2| , it is argued that 
the electron (fermion) motion becomes ballistic at a cer- 
tain quantum condition. But all fermions are known to 
be subject to scattering. It is difficult to justify the rea- 
son why the ballistic electron is not scattered by impu- 
rities and phonons, which naturally exist in nanotubes. 
Yao, Kane and Dekker [l| emphasized the importance of 
phonon scattering effects in their analysis of their data 
in Fig. [T] The Cooper pairs fli1 | in supercurrents, as is 
known, can run with no resistance (due to impurities and 
phonons) . Clearly the experiments on the currents shown 
in Fig. [T] are temperature-dependent, indicating the im- 
portance of the electron-phonon scattering effect. If the 
ballistic electron model is adopted, then the phonon scat- 
tering cannot be discussed within the model's framework. 
We must go beyond the ballistic electron model. 

If the SWNT is unrolled, then we have a graphene 
sheet, which can be superconducting at a finite tempera- 
ture. We first study the conduction behavior of graphene 
in Sec. |nl starting with the honeycomb lattice and intro- 
ducing "electrons" and "holes" based on the orthogonal 
unit cell. Phonons are generated based on the same or- 
thogonal unit cell. In Sec. Mil we treat phonons and 
phonon-exchange attraction. In Sec. 4, we construct a 
Hamiltonian suitable for the formation of the Cooper 
pairs. We derive, in Sec. 5, the linear dispersion rela- 
tion for the center-of-mass motion of the pairons. The 
pairons moving with a linear dispersion relation under- 
goes a Bose-Einstein condensation (BEC) in 2D, which is 
shown in Sec. 6. Zero-bias anomaly is discussed is Sec. 7. 

II. GRAPHENE 

Following Ashcroft and Mermin [l3[, we adopt the 
semiclassical model of electron dynamics in solids. In 



3 



(a)^. (b) 

y-axis 




FIG. 3: (a) WS unit cell, rhombus (dotted lines) for graphene. 
(b) The orthogonal unit cell, rectangle (dotted lines). 

the semiclassical (wave packet) theory, it is necessary to 
introduce a fc-vector 

k = k x e x + k y e y + k z e z , (9) 

e x , e y and e z are Cartesian orthonormal vectors since 
the fc-vectors are involved in the semiclassical equation 
of motion: 

Hk 

hk = h— = q(E + v x B), (10) 
at 

where E and B are the electric and magnetic fields, re- 
spectively, and the vector 



is the electron velocity where e is the energy. The 2D 
crystals such as graphene can also be treated similarly, 
only the z-components being dropped. The choice of the 
Cartesian axes and the unit cell is obvious for the cubic 
crystals. We must choose an orthogonal unit cell also for 
the honeycomb lattice, as shown below. 

Graphene forms a 2D honeycomb lattice. The WS unit 
cell is a rhombus shown in Fig.^a). The potential energy 
V(r) is lattice-periodic: 

V(r + R mn ) = V(r), (12) 

where 

R-mn = rna. x + na 2 (13) 

are Bravais vectors with the primitive vectors (sli,b.2) 
and integers (m, n). In the field theoretical formulation, 
the field point r is given by 

r = r' + R nm , (14) 

where r' is the point defined within the standard unit 
cell. Equation (fT2| describes the 2D lattice periodicity 
but does not establish the fc-space, which is explained 
below. 
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To see this clearly, we first consider an electron in a 
simple square (sq) lattice. The Schrodinger wave equa- 
tion is 

iftj^(r) - -^r V V( r ) + V^(r)V(r). ( 15 ) 
The Bravais vector for the sq lattice, Rmn, is 

R-mn = ma * + n&y 

— mae x + nae yi (a = lattice constant). (16) 

The system is lattice periodic: 

V sq (r + RW) =V sq (r). (17) 

If we choose a set of Cartesian coordinates (x, y) along 
the sq lattice, then the Laplacian term in Eq. (|15[) is given 
by 

V 2 ^,y)=(^ + ^)^,y). (18) 

If we choose a periodic square boundary with the side 
length Na, N — integer, then there are 2D Fourier trans- 
forms and (2D) fc-vectors. 

We now go back to the original graphene system. If we 
choose the a;- axis along either ai or a2, then the potential 
energy field V^(r) is periodic in the ^-direction, but it 
is aperiodic in the y-direction. For an infinite lattice 
the periodic boundary is the only acceptable boundary 
condition for the Fourier transformation. Then, there is 
no 2D fc-space spanned by 2D fc-vectors. If we omit the 
kinetic energy term, then we can still use Eq. (1121) and 
obtain the ground state energy (except the zero point 
energy). 

We now choose the orthogonal unit cell shown in 
Fig. EHb) . The unit has side lengths 

bi = V3a , b 2 = 3a , (19) 

where ao is the nearest neighbor distance between two 
C's. The unit cell has 4 C's. The system is lattice- 
periodic in the x- and y-directions, and hence there are 
2D fc-space. 

The "electron" ("hole") is defined as a quasi- electron 
that has an energy higher (lower) than the Fermi energy 
6f and "electrons" ("holes") are excited on the positive 
(negative) side of the Fermi surface with the convention 
that the positive normal vector at the surface points in 
the energy-increasing direction. 

The "electron" (wave packet) may move up or down 
along the y-axis to the neighboring hexagon sites pass- 
ing over one C + . The positively charged C + acts as a 
welcoming (favorable) potential valley for the negatively 
charged "electron", while the same C + acts as a hinder- 
ing potential hill for the positively charged "hole". The 
"hole", however, can move horizontally along the cc-axis 
without meeting the hindering potential hills. Thus the 



easy channel directions for the "electrons" ("holes") are 
along the y-(x-)&xes. 

Let us consider the system (graphene) at K. If we put 
an electron in the crystal, then the electron should occupy 
the center O of the Brillouin zone, where the lowest en- 
ergy lies. Additional electrons occupy points neighboring 
the center O in consideration of Pauli's exclusion princi- 
ple. The electron distribution is lattice-periodic over the 
entire crystal in accordance with the Bloch theorem. [l3j 

Carbon (C) is a quadrivalent metal. The first few 
low-lying energy bands are completely filled. The upper- 
most partially filled bands are important for the trans- 
port properties discussion. We consider such a band. The 
Fermi surface, which defines the boundary between the 
filled and unfilled fc-spaces (area) is not a circle since the 
x-y symmetry is broken. The "electron" effective mass is 
lighter in the y-direction than perpendicular to it. Hence 
the electron motion is intrinsically anisotropic. The neg- 
atively charged "electron" is near the positive ions C + 
and the "hole" is farther away from C + . Hence, the gain 
in the Coulomb interaction is greater for the "electron". 
That is, the "electron" is more easily activated. Thus, the 
"electrons" are the majority carriers at zero gate voltage. 

We may represent the activation energy difference 
by i 

ex < e 2 . (20) 

The thermally-activated (or excited) electron densities 
are given by 

n j (T) = n j e-^ k ° T , (21) 

where j = 1 and 2 denote the "electron" and "hole", re- 
spectively. The prefactor rij is the density at the high- 
temperature limit. 

III. PHONONS AND PHONON EXCHANGE 
ATTRACTION 

Phonons are bosons corresponding to the running nor- 
mal modes of the lattice vibrations. They are character- 
ized by the energy fku, where lj is the angular frequency, 
and the momentum vector ?iq, whose magnitude is 2^% 
times the wave numbers. The g-vector for phonons is 
similar to the fc-vector for the conduction electrons. The 
phonon with q represents a plane- wave proceeding in the 
q-direction. The frequency uj is connected with the q- 
vector through the dispersion relation 

oj = Lu(q). (22) 

The excitation of the phonons can be discussed based 
on the same rectangular unit cell introduced for the con- 
duction electrons. We note that phonons can be dis- 
cussed naturally based on the orthogonal unit cells. [It 
is difficult to describe phonons in the WS cell model.] 
For example, longitudinal (transverse) phonons proceed- 
ing upwards are generated by imagining a set of plates 
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each containing a number of rectangular cells executing 
small oscillations vertically (horizontally). A longitudi- 
nal wave proceeding in the crystal axis x, is represented 
by 

u q exp(— iuj q t + iq • r) = u q exp(— iuj q t + iqx), (23) 

where u q is the displacement in the x-direction. If we 
imagine a set of parallel plates containing a great number 
of ions fixed in each plate, then we have a realistic picture 
of the lattice vibration mode. The density of ions changes 
in the x-direction. Hence, the longitudinal modes are 
also called the density-wave modes. The transverse wave 
mode can also be pictured by imagining a set of parallel 
plates containing a great number of ions fixed in each 
plate executing the transverse displacements. Notice that 
this mode generates no charge-density variation. 

The Fermi velocity vf in a typical metal is of the or- 
der 10 6 ms~ 1 while the speed of sound is of the order 
10 3 ms _1 . The electrons are likely to move quickly to 
negate any electric field generated by the density varia- 
tions associated with the lattice wave. In other words, the 
electrons may follow the lattice waves instantly. Given a 
traveling normal wave mode in Eq. (|2ip , we may assume 
an electron density variation of the form: 



C q exp(— \L>j q t + iq • r). 



(24) 



Since electrons follow phonons immediately for all uj qi 
the coefficient C q can be regarded as independent of u q . 
If we further assume that the deviation is linear in the 
scalar product u q ■ q = qu^ and again in the electron 
density n(r), we then obtain 



C q = A q gw q n(r). 



(25) 



This is called the deformation potential approximation. 
[Tij The dynamic response factor Aq is necessarily com- 
plex since the traveling wave is represented by the ex- 
ponential form. Complex conjugation of Eq. (|2"4")l yields 
C*exp(iw 9 t — iq • r). Using this form we can reformu- 
late the electron's response, but the physics must be the 
same. From this consideration, we obtain 



Aq — A_ q . 



(26) 



Each normal mode corresponds to a harmonic oscilla- 
tor characterized by (q, u> q ). The displacements u q can 
be expressed as 



iq = 1 (2^) ^ ~ O 



(27) 



where (a q ,a q ) are operators satisfying the Bose 
tation rules: 



[o q , a P ] = a q a P - a p a q = <5 pq , 
[o q , a p ] = [a q , a P ] = 0. 



(28) 



Let us now construct an interaction Hamiltonian Hp, 
which has the dimensions of an energy and which is Her- 
mitian. Using Eqs. |22|) and ([23| . we obtain 



H F = [d 3 rJ2 ^[4tfu q exp(iq • r)n(r) + h.c], (29) 



where h. c. denotes the Hermitian conjugate. This Hamil- 
tonian Hp can be expressed as 



H f = 2 ( y 9 c k+q Cka i + h - c ) > 



k q 



V q = A q (h/2u, q )V 2 iq, 



(30) 



where c, are electron operators satisfying the Fermi 
anticommutation rules: 

{c k ,c k ,} = {4,4,} -0. (31) 

The Hp in Eq. ([3H)l is the Frohlich Hamiltonian. In 
the process of deriving Eq. (j30)) . we found that the Hf 
is applicable for the longitudinal phonons only. As noted 
earlier, the transverse lattice normal modes generate no 
charge density variations, making its contribution to Hp 
negligible. 

IV. THE FULL HAMILTONIAN 

Bardeen, Cooper and Schrieffer (BCS) published a his- 
toric theory of superconductivity in 1957. [16] Following 
BCS, Fujita and his collaborators developed a quantum 
statistical theory of superconductivity in a series of pa- 
pers. jT^-HII Following this theory, we construct a gen- 
eralized BCS Hamiltonian in this section. 

In the ground state there are no currents for any sys- 
tem. To describe a supercurrent, we must introduce mov- 
ing pairons, that is, pairons with finite center-of-mass 
(CM) momenta. Creation operators for "electron" (1) 
and "hole" (2) pairons are defined by 



B 



(i)t = B a)i = j j 



kitkai - C 1 C 2' 



B 



(2)f _ (2)f (2)t 
34 — c 4 c 3 • 



(32) 



We calculate the commutators among B and B\ and 
obtain 



0, 



B 



0. 



(33) 



r>U) R 0')t 
- D 12 ' -°34 



1-n^ 

„00Ji)t 

-2 L 4 



,00 



if ki = k3 and k2 = k4 
if ki = k3 and k.2 7^ k4 
if ki 7^ k3 and k2 = k4 
otherwise. 



(34) 
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Pairon operators of different types j always commute: 



V. MOVING PAIRONS 



and 



,00 



L kit L kit' 



= ifj^i, 



l 2 — c k 2 l c k 2 i 



(35) 



(36) 



represent the number operators for "electrons" (j = 1) 
and "holes" (j = 2). 

Let us now introduce the relative and net momenta 
(k, q) such that 



k = (k x - k 2 )/2 
q = ki + k 2 



k x = k + q/2 
k 2 = k + q/2. 



(37) 



Alternatively we can represent pairon annihilation oper- 
ators by 



B 



B 



'(i) 

kq 

'(2) 
kq 



B (D = (1) (1) 

- D feltfe24- — L -k+q/2i L k+q/2f 



.(2) 



C k+q/2t C -k+q/2i- 



(38) 



The prime on B will be dropped hereafter. In the 
k-q representation the commutation relations are re- 
expressed as 



^kq' ^k'q' 



= o, 



= 0, 



(39) 



The phonon exchange attraction is in action for any 
pair of electrons near the Fermi surface. In general the 
bound pair has a net momentum, and hence it moves. 
The energy w q of a moving pairon can be obtained from: 

w q a(k, q) =[e(|k + q/2|) + e(| k + q/2|)]a(k, q) 

which is Cooper's equation, Eq. (1) of his 1965 Physical 
Review Letter. The prime on the fc'-integral means 
the restriction on the integration domain arising from the 
phonon exchange attraction, see below. We note that the 
net momentum q is a constant of motion, which arises 
from the fact that the phonon exchange attraction is an 
internal force, and hence cannot change the net momen- 
tum. The pair wavefunctions a(k, q) are coupled with 
respect to the other variable k, meaning that the exact 
(or energy-eigenstate) pairon wavefunctions are superpo- 
sitions of the pair wavefunctions a(k, q). 

Equation (l4"2"| can be solved as follows. We assume 
that the energy w q is negative: 



w q < 0. 



(43) 



R 0) R Mt 

^kq I ^k'q' 



' (l ~ n-k+q/2f - n_k+q/24.) S Ji 
if k = k' and q = q' . 

M (j)t r. 

C -k+q/24. C -k'+q72i°J l 

if k + q/2 = k' + q'/2 

and - k + q/2 ^ -k' + q'/2. 



Then. 



„0')t 



(40) 



C k+q/2t C k'+q72t^' i 

if k + q/2 ^ k' + q'/2 
and - k + q/2 = -k' + q'/2. 


t otherwise. 

Using the new notations, we can write the full Hamilto- 
nian as 



k,s 



R (l)t R (2)t 
^kq ^k'q" 



o(2) R (l) 

^kq^k'q 



k' 

R (2) R (2)t 
^kq^k'q 



(41) 



This is the full Hamiltonian for the system, which can 
describe moving pairons as well as stationary pairons. 
Here, the prime on the summations indicates the restric- 
tion arising from the phonon exchange attraction, see 
below. The connection with BCS Hamiltonian [l6| will 
be discussed in Sec. IVI1 



e(|k + q/2|) + e(| - k + q/2|) - w q > 0. (44) 



Rearranging the terms in Eq. (|4"2"|) and dividing by e(|k + 
q/2|) + e(| — k + q/2|) — w q , we obtain 

a(k,q) = [ £ (|k+q/2|)+ £ (|-k+q/2|)- W9 ]- 1 C(q), (45) 
where 

I'O 



(46) 



which is fc-independent. 

Introducing Eq. (|4"5j) in Eq. (|42l) . and dropping the 
common factor C(q), we obtain 

1 =J^2 J' d3fc [^(|k + q/2|) + s(\ - k + q/2|)+ 
Hi}' 1 - (47) 

We now assume a free electron moving in 3D. The 
Fermi surface is a sphere of the radius (momentum) 



kp = (2mi£F) 



1/2 



(48) 



where mi represents the effective mass of an electron. 
The energy e(|k|) is given by 



h 2 - h 2 



(49) 
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Such a linear energy-momentum relation is valid for 
pairons moving in any dimension (D). However, the co- 
efficients slightly depend on the dimensions; in fact 

w q = wq + cq. (55) 

c/vp = 1/2 and 2/ir for 3D and 2D, respectively. 



VI. THE BOSE-EINSTEIN CONDENSATION OF 
PAIRONS 




k F +k D 



FIG. 4: The range of the integration variables (k, ( 
to a spherical shell of thickness ko- 



is limited 



The prime on the /c-integral in Eq. (|47p means the re- 
striction: 



0<e(|k + q/2|), e(| -k + q/2|) < fiwo. 



(50) 



We may choose the polar axis along q as shown in Fig. |U 
The integration with respect to the azimuthal angle sim- 
ply yields the factor 2tt. The fc-integral can then be ex- 
pressed by 

(n~%\3 r 71 / 2 r k F +k D -±q cos 6 

^L=4n ddsm9 



kp-\- hq cos t 



k 2 dk 



\w q \ + 2e k + (4m 1 )- 1 q 2 



mihuJuk v 1 . 



(51) 



(52) 



After performing the integration and taking the small-g 
and small- (&d/A;f) limits, we obtain 



w q =w a + —q, 



(53) 



where the pairon ground-state energy wq is given by 



w 



exp[2/uoAf(0)] - 1 ' 



(54) 



As expected, the zero-momentum pairon has the lowest 
energy wq. The excitation energy is continuous with no 
energy gap. Equation (|53p was first obtained by Cooper 
(unpublished), and it is recorded in Schrieffer's book [22| . 
Eq. (2.15). The energy w q increases linearly with mo- 
mentum (magnitude) q for small q. This behavior arises 
from the fact that the density of states is strongly re- 
duced with the increasing momentum q and dominates 
the q 2 increase of the kinetic energy. The linear disper- 
sion relation means that a pairon moves likes a massless 
particle with a common speed v-p/2. This relation plays 
a vital role in the B-E condensation of pairons (see next 
section). 



In Sec. IIV1 we saw that the pair operators (B, B<) ap- 
pearing in the full Hamiltonian H in Eq. (|4T|) satisfy 
rather complicated commutator relations in Eqs. (|39p and 
(HO)) . In particular part of Eq. 



By 



kO 



,t J 
'-kt c kj. 







(56) 



reflect the fermionic natures of the constituting electrons. 
Here, B^ = 6^ represents creation operator for zero 
momentum pairons. BCS [16j studied the ground-state 
of a superconductor, starting with the reduced Hamilto- 
nian Hq, which is obtained from the Hamiltonian H in 
Eq. (|41[) by retaining the zero momentum pairons with 



q = 0, written in terms of b by letting B 



00 

kO 



dp) 



-(l) 6 (l)t 6 (l) 
k °k "k 



b$b£) +b w b w 



k °k' "k 



(57) 



Here, we expressed the "electron" and "hole" kinetic en- 
ergies in terms of pairon operators. The reduced Hamil- 
tonian Hq is bilinear in pairon operators (6, w), and can 
be diagonalized exactly. BCS obtained the ground-state 
energy Eq as 



E = hu) D Af(0)w o , 



(58) 



where Af(0) is the density of states at the Fermi energy. 
The wq is the ground-state energy of the pairon, see 
Eq. (fS~4")) . Equation ([55)1 means simply that the ground 
state energy equals the numbers of pairons times the 
ground-state energy wo of the pairon. Our Hamiltonian 
H in Eq. (|¥T1l is reduced to the original BCS Hamiltonian 
(see Ref. [16(, Eq. (24)). There is an important differ- 
ence in the definition of "electron" and "hole" here. BCS 
called the quasi-electron whose energy is higher (lower) 
than the Fermi energy ep, the "electron" ("hole"). In our 
theory the "electrons" ("holes") are defined as quasiparti- 
cles generated above (below) the Fermi energy and cir- 
culates counterclockwise (clockwise) viewed from the tip 
of an external magnetic field vector B. They are gener- 
ated, depending on the energy contour curvature signs. 
For example, only "electrons" ("holes") are generated for 



s 



a circular Fermi surface with negative (positive) curva- 
ture whose inside (outside) is filled with electrons. Since 
the phonon has no charge, the phonon exchange cannot 
change the net charge. The pairing interaction terms in 
Eq. (|4Tj) conserve the charge. The term — VqB^Ib^S s , 

where vq = |VqV^| {TiloqA^ 1 , A — sample area, gener- 
ates a transition in the "electron" states. Similarly, the 
exchange of a phonon generates a transition in the "hole" 

(2H (2)t 

states represented by — vqB^^ B^, s . The phonon ex- 
change can also pair-create or pair-annihilate "electron" 
("hole") pairons, and the effects of these processes are rep- 
resented by -VoB^B^, ^B^B^, as shown in 
Feynman diagrams in Figs. 0(a) and0(b). At K the sys- 
tem must have equal numbers of — (+) zero-momentum 
(ground) pairons. 

To describe a supercurrent, we must introduce moving 
pairons. We now show that the center-of-masses of the 
pairons move as bosons. That is, the number operator of 
pairons having net momentum q 



kq 



(59) 



have the eigenvalues 

n q = 0,l,2,.... (60) 

The number operator for the pairons in the state (k, q) 
is 

"kq = ^ q S kq = 4+ q /2 C -k+q/2 C -k+q/2Ck+ q /2' ( 61 ) 

where we omitted the spin indices. Its eigenvalues are 
limited to zero or one: 



'kq 



= or 1. 



(62) 



To explicitly see this property in Eq. (|60p . we introduce 



and obtain 



[^q,»q]=E( 1 
k 

[n q ,B q ] =B\. 



- «k+Iq - n -k+iq ] B kq = #q, 



(64) 



Although the occupation number n q is not connected 
with Bq as n q ^ _B q _B q , the eigenvalues n q of n q satis- 
fying Eq. (|64|) can be shown straightforwardly to yield 



n q = 0,l,2,.. 



with the eigenstates 

|0>, |l)=Bt| ), |2)=flt B t| 



(65) 



(66) 



In summary, pairons with both k and q specified are 
subject to the Pauli exclusion principle, see Eq. (|()2"j) . Yet, 



the occupation numbers n q of pairons having a CM mo- 
mentum q are 0, 1, 2, ... . 

The most important signature of many bosons is the 
Bose-Einstein Condensation (BEC). Earlier we showed 
that the pairon moves in 2D with the linear dispersion 
relation, see ([55)1 : 



u p = w + (2/tt)v f p = w Q + cp, 



(67) 



where we designated the pairon net momentum (magni- 
tude) by the more familiar p rather than q. 

Let us consider a 2D system of free bosons having a 
linear dispersion relation: e = cp, c = (2/tt)vf. The 
number of bosons, N, and the Bose distribution function 



/ B (e;/3,M) 



= / B ( £ )(>0) (68) 



are related by 

n = J2Mc p ;P^) = n + J2'M e p)> ( 69 ) 

P p 

£p>0 

where ix is the chemical potential, j3 = (fee? 1 )" 1 , and 

N = (e-^ 1 - l)- 1 (70) 

is the number of zero-momentum bosons. The prime on 
the summation in Eq. (|69p indicates the omission of the 
zero-momentum state. For notational convenience, we 
write 



e = cp — (2/tt)vfp (> 0). 



(71) 



We divide Eq. (|69| by the normalization area L 2 , and 
take the bulk limit: 

N -> oo, L -> oo while NL~ 2 = n. (72) 

We then obtain 

n-n = J rf2 P/B(e), (73) 

where no = No/L 2 is the number density of zero- 
momentum bosons and n the total boson density. After 
performing the angular integration and changing integra- 
tion variables, we obtain from Eq. (|73ll : 



2<Kh 2 c 2 l3 2 (n-no) 



dx- 



o A" 
(x = fie), 



where the fugacity 



X = e^ (< 1). 



(74) 



(75) 



is less than unity for all temperatures. After expanding 
the integrand in Eq. (|74[) in powers of Xe~ x (< 1), and 
carrying out the x-integration, we obtain 



n x = n-n 



fcjT 2 2 (A) 
2ttTi 2 c 2 ' 



(76) 
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FIG. 5: Feynman diagrams representing (a) pair-creation of ± ground pairons from the physical vacuum, and (b) pair annihi- 
lation. The time is measured upwards. 



fc=i 



A- 2 



(0 < A < 1) 



(77) 



Equation (|76[) gives a relation among A, n, and T. 

The function 02(A) monotonically increases from zero 
to the maximum value 



2 (1) = 1-645 



(78) 



as A is raised from zero to one. In the low-temperature 
limit, A = 1, </>2(A) = 02(1) = 1-645, and the den- 
sity of excited bosons, n x , varies as T 2 as seen from 
Eq. Q76p. This temperature behavior of n x persists 
as long as the right-hand-side (r.h.s.) of Eq. ([75)) is 
smaller than n; the critical temperature T c occurs at 



fc B T 2 2 (l)/2^ 



2 „2 



c . Solving this, we obtain 



fc B T r = 1.954am 1 / 2 ( = 1.247i?W /2 



(79) 



The BEC of pairons moving in 2D occurs at a finite 
temperature. This appears to contradict with Hohen- 
berg's theorem (no long range order in 2D). But this 
theorem is proved under the assumption of the f-sum 
rule arising from the mass conservation. The pairons 
move massless with the linear dispersion relation [see Eq. 
(|71])]. and hence they are not subject to Hohenberg's the- 
orem, [m 

If the temperature is raised beyond T c , the density of 
zero momentum bosons, uq, becomes vanishingly small, 
and the fugacity A can be determined from 



fc B r 2 2 (A) 

2irh 2 c 2 



T>T C 



(80) 



In summary, the fugacity A is equal to unity in the 
condensed region: T < T c , and it becomes smaller than 
unity for T > T c , where its value is determined from 
Eq. ([HDD. 

Formula (|79l) for the critical temperature T c is distinct 
from the famous BCS formula 



3.53A:bT c = 2A , 



(81) 



where A is the zero temperature electron energy gap 
in the weak coupling limit. The electron energy gap 



A(T) and the pairon ground-state energy wq both de- 
pend on the phonon-exchange coupling energy parame- 
ter vq, which appears in the starting Hamiltonian H in 
Eq. (|41[). The energy wq is negative (bound-state en- 
ergy). Hence, this wq cannot be obtained by the pertur- 
bation theory. The connection between wq and Do is very 
complicated. This makes it difficult to discuss the critical 
temperature T c based on the BCS relation ([ST]) . Unlike 
the BCS formula, formula (f79")) is directly connected with 
the measurable quantities: the pairon density Uq and the 
Fermi speed v-p. 

We emphasize here that both formulas (|79l) and (|8Tj) 
were derived, starting with the Hamiltonian H in Eq. (|lTj) 
and follo wing statistical mechanical calculations, see the 
reference [26] for details. 



VII. ZERO BIAS ANOMALY 

The unusual current-dips at zero bias in Fig. [1] is often 
called the zero bias anomaly (ZBA). This effect is clearly 
seen in (a) low resistance contacts LRC sample. The 
differential conductance dl/dV increases with increasing 
bias, reaching a maximum at V ~ 100 mV. With a further 
bias increase, dl/dV drops dramatically. See (a), the 
upper panel in Fig. [TJ We will show that the ZBA arises 
from the break-down of the superconducting state of the 
system. 

With no bias, the nanotube's wall below ~ 100 K is 
in a superconducting state. If a small bias is applied, 
then the system is charged, positively or negatively de- 
pending on the polarity of the external bias. The applied 
bias field will not affect the neutral supercurrent but can 
accelerate the charges at the outer side of the carbon 
wall. The resulting normal currents carried by conduc- 
tion electrons are scattered by impurities and phonons. 
The phonon population changes with temperature, and 
hence the phonon scattering is temperature-dependent. 
The normal electric currents along the tube length gen- 
erate circulating magnetic fields, which destroy eventu- 
ally the supercurrent running in the wall at a high enough 
bias. Thus, the current I (/J.A) versus the voltage V (mV) 
is non-linear near the origin because of the supercurrents 
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running in the wall. The differential conductance dl/dV 
is very small and nearly constant (superconducting) for 
V < 10 mV in the HRC sample, see the lower panel in 
Fig. [1] We stress that if the ballistic electron model Q is 
adopted, then the scatterings by phonons cannot be dis- 
cussed. The non-linear IV curves below 150 K mean that 
the carbon wall is superconducting. Thus, the clearly 



visible temperature effects for both LRC and HRC sam- 
ples arise from the phonon scattering. We assumed that 
the system is superconducting below ~ 150 K. The ZBA 
arises only from the superconducting state. The super- 
conducting critical temperature T c must then be higher 
than 150 K. An experimental check of T c is highly desir- 
able. 
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